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Abstract 

p-Aminobenzoic acid (pABA) plays important roles in a wide variety of metabolic processes. Herein we report the synthesis, 
theoretical calculations, crystallographic investigation, and in vitro determination of the biological activity and phytotoxicity 
of the pABA salt, 2-hydroxyethylammonium p-aminobenzoate (HEA-pABA). The ability of neutral and anionic forms of pABA 
to interact with TIR1 pocket was investigated by calculation of molecular electrostatic potential maps on the accessible 
surface area, docking experiments, Molecular Dynamics and Quantum Mechanics/Molecular Mechanics calculations. The 
docking study of the folate precursor pABA, its anionic form and natural auxin (indole-3-acetic acid, IAA) with the auxin 
receptor TIR1 revealed a similar binding mode in the active site. The phytotoxic evaluation of HEA-pABA, pABA and 2- 
hydroxyethylamine (HEA) was performed on the model plant Arabidopsis thaliana ecotype Col 0 at five different 
concentrations. HEA-pABA and pABA acted as potential auxin-like regulators of root development in Arabidopsis thaliana 
(0.1 and 0.2 mM) and displayed an agravitropic root response at high concentration (2 mM). This study suggests that HEA- 
pABA and pABA might be considered as potential new regulators of plant growth. 
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Introduction 

The interest in carboxylate compounds in general and benzoic 
acid derivatives in particular continues to grow because of their 
chemical and biological properties. Benzoic acids and their 
derivatives are important structural elements for many natural 
products, being involved in various physiological processes in 
plants. Among these compounds, p-aminobenzoic acid (pABA) has 
displayed biological, medicinal and industrial interest. ^ABA is a 
compound of high biological significance, it is present in plant and 
animal tissues and is sometimes referred to as bacterial vitamin H t , 
B x , or Bio [1]. J^ABA is a well-known precursor of folic acid [2], 
and has recently been identified as a coenzyme Q_ precursor [3,4], 
which opens new possibilities for its application in medicine. It is a 
building block used in design of drugs and frequendy found as a 
structure moiety in drugs [5]. pABA exhibits a wide range of 
therapeutic uses as antioxidant [6,7], antibacterial [8,9], antimu- 
tagenic [10], anticoagulant [11,12], fibrinolytic and immunomod- 
ulating agent [13], protective drug against UV-irradiation [14-16] 
and in diagnostic tests for the state of the gastrointestinal tract 
[17,18]. 

Recent studies in agriculture report the role of pABA in plant 
thermotolerance [19] and its action as chemical inducer of 
systemic acquired resistance against plant pathogens [20]. 
Considered to be in the B-complex vitamin family and charac- 



terized by safety and cost effectiveness, ^ABA is frequendy added 
along with plant hormones in nutriment media as a stimulator of 
seed germination [21]. To the best of our knowledge, the role of 
folate precursor pABA as a plant growth modulator has not been 
investigated so far. 

Besides the plethora of properties noted above, ^ABA is a 
versatile reagent for structure extension through linear hydrogen 
bonding associations, through both the carboxylic and amine 
functional groups [22]. Keeping in view the structural and 
biological diversity of pABA and in connection with our interest 
in supramolecular chemistry of carboxylate salts [23,24], here we 
report the synthesis, theoretical calculations, crystallographic 
investigation, and in vitro phytotoxic activity of 2-hydroxyethy- 
lammonium jf?-aminobenzoate (HEA-^ABA). As a base, we chose 
the natural amino alcohol 2-hydroxyethylamine (HEA), a suitable 
model of alkanolamine. This compound is cheap, commercially 
available, environmentally tolerable and has hydrogen-bond 
donor sites. On the other hand, HEA is an essential component 
of cell membranes and an important metabolite in plants for the 
synthesis of choline and membrane lipids, such as phosphatidyl- 
ethanolamine and phosphatidylcholine [25] . 

The preparation of multicomponent organic crystals (solvates, 
hydrates, co-crystals, and salts) from such multifunctional mole- 
cules as benzoic acids and alkanolamine display interesting 
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networks and a wide range of different biological properties 
[23,24,26-28]. The possible binding modes of pABA, its anionic 
form and the natural auxin IAA bound to the auxin receptor TIR1 
have been evaluated by docking experiments and refined by 
Molecular Dynamics (MD) and Quantum Mechanics/Molecular 
Mechanics (QM/MM) calculations. The present work continues 
our recent investigations of new structural aspects of alkanolamine 
salts [24], reports the full details of the X-ray crystallographic 
structure of HEA-/ABA and introduces this compound and 
corresponding acid as new potential regulators of plant growth. 

Materials and Methods 

General 

/ABA and HEA used for the synthesis of YLEA-pABA were of 
analytical grade and purchased from Fluka AG (Buchs SG). HEA 
was freshly distilled before any use. FT-IR spectra in the range 
4000-400 cm" 1 were recorded with a JASCO - FT/IR-4200 
spectrometer, using the KBr pellet technique with a resolution of 
4.0 cm 1 and a scanning speed of 16 mm s . The optical 
properties of HEA-/?ABA were examined by using a UV— Vis 
spectrophotometer at room temperature. UV— visible spectra were 
recorded in the 190-800 nm range with PERKIN-ELMER 
LAMBDA 12 UV— Vis spectrometer. Melting point was found for 
finely purifying compound (by repeated recrystallization) using a 
Boetius instrument. 

Synthesis and characterization of HEA-pABA 

HEA-/) ABA was prepared by the reaction of HEA andp ABA in 
a 1:1 molar ratio. Under vigorous stirring, freshly distilled HEA 
was added drop wise to the solution of the acid in acetone, at 
reflux for 2 h to complete the reaction. The formed salt was 
precipitated in a white crystalline state and, after cooling at room 
temperature, collected by filtration, washed with cold acetone, and 
dried under vacuum for 3 h. This reaction had a 90% yield. The 
salt was recrystallized to obtain suitable crystals for X-ray analysis. 
Crystals were grown by slow evaporation of ethanol solution at 
room temperature. The purity of synthesized salts was 99.3% 
(±0.8%), established by spectrophotometric method. This con- 
sisted in determining the specific extinction (absorbance at X max of 
a solution 1 gl in a cell of 1 cm width) of the free acid (s a ) and 
the corresponding salt (s s ), for the UV maxima: 

s s -M s 
purity = ; - -100, 

s a -M a 

where M„ and M a are the molecular weights for the salt and free 
acid, respectively. The elemental analysis was in agreement with 
the expected stoichiometry. The FT-IR spectrum (KBr pellets 
prepared by grinding a 5 — 10 mg sample with 100 mg KBr) was 
consistent with salt formation. 

C 9 H 14 N 2 0 3 , (198.22), m.p. 147-150°C; X max = 265.19 nm; FT- 
IR spectra (KBr pelet, cm" 1 ): 3392,3227,2132, 1591, 1540, 1384, 
1361, 1073, 1016; calcd. (%): C 54.53, H 7.12, N 14.13; found 
(%): C 54.51, H 7.03, N 14.10. 

X-ray measurement and refinement 

X-ray data for HEA-/>ABA were collected at room temperature 
on a Siemens P3/PC diffractometer (CuK a -radiation). The 
crystallographic data and the experimental details are summarized 
in Table 1. Structure solution and refinement were carried out 
using the SHELX-97 program [29]. Non-hydrogen atoms were 
refined anisotropically. The hydrogen atoms were placed in 
calculated positions with their isotropic displacement parameters 



riding on those of parent atoms and the H-atoms of NH 2 group of 
HEA-/ABA were found from differential Fourier maps. The 
geometric parameters for H-bonds are listed in Table 2. The 
geometric parameters were calculated and the figures were drawn 
with the use of the PLATON program [30] . Bond lengths (A) and 
angles (°) for HEA-/ABA are listed in Table 3. The hydrogen 
atoms that are not involved in the hydrogen bonding were omitted 
from the representation of the crystal packing. Crystallographic 
data for HEA-/>ABA were deposited in the Cambridge Crystal- 
lographic Data Center (CCDC 936586). 

Docking and QM/MM calculations 

The natural auxin (indole-3-acetic acid, IAA) and /ABA, the 
last one both in the neutral and anionic form, were docked in the 
auxin binding site as identified in the crystal structure of the TIR 1 
protein-auxin complex retrieved from protein data bank (PDB id 
2P1P) [31] using AutoDock 4.2 software [32]. This software 
considers only the polar hydrogen atoms during the docking 
process. A Lamarkian genetic algorithm [33] was employed for the 
docking simulation, performing 100 independent runs per 
molecule. In each run, a population of 50 individuals evolved 
along 27000 generations and a maximum number of 25 million 
energy evaluations were performed. The best fit (lowest docked 
energy) solutions of the 100 independent runs were stored for 
subsequent analysis. 

Docking results were refined by performing molecular dynamics 
simulations using the AMBER 11 program suite [34], the 
AMBER99SB force field [35] for the protein and the GAFF force 
field for ligands [36]. Each of the three protein-ligand complexes 
was first submitted to a geometry optimization consisting of 1000 
steps of steepest descent, followed by 9000 steps of conjugate 
gradient, energy minimization. The systems were then equilibrat- 
ed for 100 ps in the NVT ensemble and for additional 100 ps in 
the NPT ensemble. Ten ns long molecular dynamics simulations 
were then run in the NPT ensemble at 300 K and 1 atm. Water 
solvent was modelled with the TIP4P water model [37]. A cutoff of 
12 A was used for non-bonded VdW interactions, while electro- 
static interactions were treated with the Particle Mesh Ewald 
approach [38]. 

QM/MM calculations were performed on the final structures 
obtained from MD simulations, using the two-layer ONIOM 
model [39-41] as implemented in Gaussian09 [42]. Owing to the 
high complexity of the examined system, QM/MM calculations 
were performed in vacuo, including only the ligands and the 
residues that showed the most favorable interactions with the 
ligands during the MD simulations (namely Arg403, Ser438 and 
Leu439) in the QM layer. The protein was included in the MM 
layer and its geometry was kept fixed. The QM part was treated at 
Density Functional Theory (DFT) level with the B3LYP functional 
[43,44] and the 6-3 1+G* basis set. Protein atom charges were 
assigned using the Amber Force Field [34]. The mutual 
electrostatic interaction between the ligand and the protein was 
treated using the electrostatic embedding approach [41]. In this 
model, the point charges of the MM layer are included in the 
Hamiltonian of the QM layer, introducing a polarization of the 
wave function. 

Biological activity 

Plant material and seed germination. Arabidopsis thali- 
ana, ecotype Columbia 0 (Col 0) seeds were surface sterilized by 
treatment with a 5% (w/v) calcium hypochlorite solution and 
0.02% (v/v) Triton X-100 in 80% ethanol (EtOH), for 10-12 
minutes at 25-28°C, by continuous shaking. They were then 
rinsed twice with 80% EtOH, washed with 100% EtOH and 
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Table 1. Crystal Data and Structure Refinement for HEA-pABA. 





Compound 


HEA-pABA 


Empirical formula 


CgH, 4 N 2 0 3 


Formula weight 


198.22 


Wavelength 


1.54184 


Crystal system 


Orthorhombic 


Space group 


Pbca 


a, A 


1 1 .898(2) 


fa, A 


8.3140(17) 


C A 


20.790(4) 


a, deg 


90 


P, deg 


90 


y, deg 


90 


V, A 3 


2056.5(7) 


Z 8 


Density (calc), mg/m 3 


1.280 


Absorption coefficient, mm-1 


0.808 


F(000) 


848 


Crystal size, mm3 


0.14x0.10x0.08 


0 range for data collection, deg 


4.25-70.06 


Index ranges 


-14S/1S14, -lOsfel, -25</<24 


Reflections collected/unique 


1952/1952, fi(int) = 0.0000 


Completeness (%) 


99.9 (6 = 70.06) 


Data/parameters 


1952/135 


GOF on F2 


0.998 


Final R indices (/>2a (/)) 


fil =0.0467, wR2 = 0.0930 


ft indices (all data) 


fil =0.1101, w/?2 = 0.1105 


Largest diff. peak and hole, e.A~ 3 


0.196/-0.193 


doi:1 0.1 371 /journal.pone.01 01 892.W01 



finally with sterile distilled water [45]. Thirty sterile seeds were 
sown on every square Petri plate (12x12 cm) contained: 80 ml of 
sterile solid MS agar medium (modified from [46]), HEA-pABA, 
pABA or HEA at five concentrations (0.02, 0.1, 0.2, 1 and 2 mM) 
and IAA at five concentrations (0.01, 0,1, 1, 10 and 100 nM). The 
salt and IAA were dissolved in 5 mM MES (2-pV'-morpholino] 
ethanesulfonic acid) buffer, while pABA, which is slightly soluble in 
water, was dissolved in MES buffer - EtOH solution, with a final 
EtOH concentration of 0.2% (v/v). Our previous research on 
Arabidopsis thaliana (Col 0) [28] and Cucumis sativus L. [27] 



seeds show that 0.2% (v/v) ethanol concentration in culture media 
does not alter germination rate of seeds, but inhibit slowly primary 
root length and lateral root formation, as it is confirmed by some 
studies reported in the literature [47,48]. For clarity, control 
corresponding to pABA was noted CONTROL 0.2% EtOH and 
for HEA-pABA and HEA, CONTROL. HEA-pABA, pABA or 
HEA solved in the solvents mentioned above were sterilized by 
filtration (Millipore) and added to molten control agar (~50°C). 
Five repetitions for each compound concentration were per- 
formed. The plates were stratified horizontally in darkness at 4°C 



Table 2. Hydrogen-bonding geometry (A) for compound HEA-pABA. 





No. of HB 


D-H-A 


d(D H), A 


d(H A), A 


d(D A), A 


Z(DHA), deg. 


Symmetry transformation for H-acceptor 


1 


N(1)-H -0(3) 


0.89 


1.91 


2.801(3) 


175 


-x+2, -y, -z+2 


2 


N(1)-H~N(2) 


0.89 


2.22 


3.051(4) 


155 


x, -y+1/2, z+1/2 


3 


N(1)-H -0(2) 


0.89 


1.87 


2.752(3) 


171 


x, y, z 


4 


0(1)-H "0(3) 


0.82 


1.90 


2.715(3) 


176 


x, y, z 


5 


N(2)-H -0(1) 


0.91(3) 


2.29(3) 


3.087(3) 


146 


-x+2, y-1/2, -z+1/2 


6 


N(2)-H -0(3) 


0.82(3) 


2.20(3) 


3.025(3) 


174 


x+1/2, y, -z+1/2 


7 


C(2)-H 0(2) 


0.97 


2.41 


3.247(3) 


145 


x-1/2, -y+1/2, -z+1 



doi:1 0.1 371 /journal.pone.01 01 892.t002 
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Table 3. Bond lengths (A) and angl 


ss (°) for HEA-pABA. 




N(1)-C(1) 


1.481(3) 


0(1)-C(2) 


1 .409(3) 


C(1)-C(2) 


1.495(4) 


N(2)-C(6) 


1 .404(3) 


0(2)-C(9) 


1.257(3) 


0(3)-C(9) 


1.281(3) 


C(3)-C(4) 


1.384(3) 


C(3)-C(8) 


1.391(3) 


C(3)-C(9) 


1.492(3) 


C(4)-C(5) 


1.377(3) 


C(5)-C(6) 


1.387(4) 


C(6)-C(7) 


1 .389(4) 


C(7)-C(8) 


1.379(3) 


N(1)-C(1)-C(2) 


112.1(3) 


0(1)-C(2)-C(1) 


113.9(3) 


C(4)-C(3)-C(8) 


117.8(2) 


C(4)-C(3)-C(9) 


121.8(2) 


C(8)-C(3)-C(9) 

^-y j / ^-v-v 


120.3(3) 


C(5)-C(4)-C(3) 


121.4(3) 


C(4)-C(5)-C(6) 


120.7(3) 


C(5)-C(6)-C(7) 


118.3(2) 


C(5)-C(6)-N(2) 


121.5(3) 


C(7)-C(6)-N(2) 


120.1(3) 


C(8)-C(7)-C(6) 


120.7(3) 


C(7)-C(8)-C(3) 


121.1(3) 


0(2)-C(9)-0(3) 


122.9(3) 


0(2)-C(9)-C(3) 


118.6(3) 


0(3)-C(9)-C(3) 


118.5(3) 


doi:1 0.1 371 /journal.pone.01 01 892.W03 



for 2 days to break dormancy and synchronize germination, and 
then grown in a climate chamber at 22°C in a vertical position 
under a photoperiod of 16 h light (220 |J,mol m s ) and 8 h 
darkness. 

Data analysis 

After 5 and 10 days, respectively, images of the Petri dishes were 
generated by scanning the individual plate on an EPSON 
Perfection 4870 Photo scanner at 600 DPI. Obtained images 
were digitized using the ImageJ software package, and primary 
root length was measured for each seedling. After 10 days, plants 
were stained with 0.1% (w/v) safranin for 1-2 h, then rinsed with 
distilled water and 100, 80, 50, 30, 15% EtOH, respectively, and 
mounted in 50% (v/v) glycerol. All lateral roots emerging from the 
primary one were observed using a Nikon YS 1 00 microscope 1 0 x 
objective, and were taken into account for lateral root number 
data. 

Statistical analysis 

For all experiments, the overall data were statistically analysed 
by SYSTAT 10.0 program. All values are expressed as a mean ± 
S.E. One-way ANOVA analyses, followed by i-test using the 
Bonferroni correction, were used to compare the data obtained for 
primary root length, lateral root number at 5 and 10 days for 



treatments (HEA-^ABA, pABA and HEA) at five different 
concentrations and controls. Different letters are used to indicate 
means that differ significantly (P<0.05). 

Results and Discussion 

Synthesis and characterization of HEA-pABA 

HEA-/fABA is a stable solid salt in accordance with the 
appropriate "rule of three". The value of ApKa (the difference 
between the pKa value for HEA and pAB A) was estimated to be > 
3, that determines the extent of proton transfer. UV spectroscopic 
data obtained in 0. 1 M NaOH showed similar X max values for 
both salt and acid, that confirms the existence of the same anion in 
both compounds (see Figure SI). The cut off wavelength for HEA- 
pABA was found to be 265.19 nm, corresponding to the gap 
energy of 4.68 eV, which is typical of insulating materials. There is 
no absorption in the entire visible region. This is the most 
favourable characteristic for a non-linear optical material. 

FT-IR analysis indicates salt formation by the presence of bands 
arising from the asymmetric and symmetric vibrations of the 
COO - group occurring at 1591 cm - and 1384 cm - , which do 
not exist in the spectra of free acid and by the absence of vC — O 
and vC-OH bands at 1673 cm 1 and 1291 cm , characteristic 
for -COOH group [49,50]. At the same time, the band of the 
deformation vibrations of the N-H bonds is observed at 
1540 cm - . A supplementary proof of the salt formation was the 
appearance of C — O vibrations between 1100 and 1000 cm , 
belonging to — CH 2 OH from HEA and a broad band near 
2100 cm -1 , which often appears in amine salt [50] (see Figure S2). 
The -NH 3 + vibrations appears at 3200-2800 cm" 1 and 5NH 3 + at 
2400-2800 cm . The changes of wavenumbers of the bands of 
carboxylate group in the case of salt were discussed comparing to 
the acid. From the spectroscopic data and the elemental analysis it 
is concluded that pABA reacted with HEA with a 1 : 1 ratio. 

Crystal structure analysis of HEA-pABA 

The salt unit of HEA-^ABA serves as building block for the 
supramolecular architecture and in the crystal it is self-assembled 
via ionic N-H "O and normal O-H "O hydrogen bonds (HB) 
(Figure 1, Table 2). The anion in HEA-^ABA forms a practically 
planar system because the dihedral angle between the least-square 
plane of the phenyl ring^4(C3C4C5C6C7C8) and the least-square 
plane of the COO - group is equal to 5.8°. The cation adopts the - 
Syn-Clinal conformation, [30] the N1C1C201 torsion angle in 
HEA-jf)ABA is equal to —68.8°. In the studied compound, the 
nitrogen atom in amino group has practically a pyramidal 
configuration because its valence angles are equal to 111 .4, 
1 14.9 and 1 17.2°. It means that the non-bonding electron pair of 
N2 could be considered as a fourth substituent on a sp 3 hybridized 
nitrogen atom. Thus, the pyramidal geometry of the amino group 
allows forming three hydrogen bonds in the crystal. Indeed, the 
building blocks are assembled into 2-D layers parallel to (100) 
plane (Figure 2 A, Table 2) via Nl "03 (1) HB and two hydrogen 
bonds with participation of amino group Nl "N2 (2) and N2 " Ol 
(5). These layers are further related along a axis by a glide plane 
perpendicular to [00 1] with the glide component [0.5 0 0] and a 
twofold screw axis with direction [100] and screw component [0.5 
0 0] due to N2 03 (6) and C2 02 (7) HB respectively (Figure 2B, 
Table 2). According to the criterion proposed by Spek [30], to 
establish the presence of Jt— 71 stacking interactions between phenyl 
rings Jl and J2 (Cg(Jl)---Cgg2)<6.0 A, P<60.0°, where Cg are 
the respective centroids of the rings and P is the angle between the 
Cg(Jl)Cg(J2) vector and the normal to the Jl ring), there is a 
stacking interaction between the phenyl rings A and A (— x, 1/2+y, 
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1/2— z) within the layers of HEA-jf)ABA. The distance between 
their centroids is in fact equal to 5.062 A and (3 = 16.7°. There is 
an effect of hydrogen bonding on the geometry of carboxylate 
groups. While the carboxylate group in HEA-^)ABA participates as 
hydrogen-bonding acceptor, the C-O bond lengths vary signifi- 
cantly with the number and type of hydrogen-bonding donors 
linked to the oxygen atom. In the absence of hydrogen bonding 
and other electronic perturbations, the C-O bond lengths should 
be equal because of electron derealization [23,24]. Formation of 
single or multiple hydrogen bonds at one oxygen atom should 
cause the associated C-O bond to lengthen. In the studied 
compound, the difference in C9-02 (1.257(3) A) and C9-03 
(1.281(3) A) bond lengths in carboxylate group is equal to 0.024 A 
(Table 3). The transfer of proton in HEA-pABA has led to a 
significant polarization of the anion in comparison with pABA 
because the calculated dipole moments using HF method [5 1] for 
these systems are equal to 13.952 and 3.13 Debye, respectively, 
and dipole vectors in both molecules lie very close to C3-C9 bond. 

Molecular properties of HEA-pABA and pABA compared 
to classical auxin molecules 

Beside the importance of salt unit of HEA-pABA as building 
block for supramolecular architecture, we found a new role of this 
compound and its precursor (pABA) on plant development as 
potential regulators of plant growth with an auxin-like effect. 
Auxins are defined by a set of physiological actions: cell division, 
elongation, formation of lateral and adventitious roots and 
gravitropism [52-54] . A large number of compounds with auxin 
and anti-auxin activity were identified so far as a result of 
structure-activity relationship investigations. More than 200 
compounds with auxinic action and with different chemical 
structures have two common features critical for auxin activity: a 
planar aromatic ring structure and a carboxyl group-containing 
side chain [55]. The activity of both natural and synthetic auxins 
depends strongly on side groups binding to the ring structure. 
Through such substitutions, it is suggested that a positive charged 
area at specific distance from the carboxylic group is generated. 
The distance between the negative charge on the carboxyl group 
and the partial positive charge on the ring has been claimed to be 
critically important for detectable activity. According to some 
authors, the optimal distance (OD) appears to be about 5.0-5.5 A 
[56,57]. The ODs in HEA-pABA andpABA were compared with 
those of classical auxin molecules: indole-3-acetic acid (IAA), 2,4- 
dichlorophenoxyacetic acid (2,4-D) and 1 -naphthaleneacetic acid 
(1-NAA), where IAA is the most important member of the auxin 
family and the natural one, the others being synthetic auxin 
analogs. We performed a search through the Cambridge 
Structural Database (version 5.34) [58,59] and found that the 
crystal structures of these compounds indeed have an OD ranging 
from 4 to 6 A. We observed 1 1, 4, 7 hits for IAA, 1-NAA and 2,4- 
D where the ODs range from 4.3 to 5.5, 4.3 to 4.6 and 5.1 to 
5.3 A, respectively. The presence of several conformers for IAA, 1- 
NAA and 2,4-D occurs due to rotation of molecules' moieties on 
bonds Cl(0)-C2 and C2-C3 (see Figure S3). However, only for 
the conformers of 2,4-D family the ODs satisfy to specific distance 
between the negative charge on the carboxyl group and the partial 
positive charge on the ring. Therefore, this result brought us to 
look for the low energy conformer for IAA and 1 -NAA, and then 
to check the above mentioned criteria for OD values. This was 
achieved by using the module Conformational Search implement- 
ed in HyperChem 6.03 [60]. The method involves the use of the 
directed scheme as introduced in the Monte Carlo multiple 
minimum approach [61]. This method seeks to uniformly sample 
low-energy regions by cycling through all previously accepted 




Figure 1. ORTEP drawing for HEA-pABA. The atomic labeling and 
charge-assisted hydrogen bonds are represented by dashed lines. 
Thermal ellipsoids are shown with the 50% probability level. 
doi:1 0.1 371/journal.pone.01 01 892.g001 

conformations (from lowest energy to highest, in order) when 
selecting each initial structure. In this process, the initial structures 
from series conformers of LAA and 1-NAA were modified by 
variation of the torsion angles along the C1-C2 and C2-C3 bonds 
followed by energy minimization of these angles using the MM+ 
force field. Then the point charges of atoms were calculated using 
HF method [51] for found low calculated energy conformers of 
IAA, 1-NAA and crystal structures of HEA-pAB A, p ABA [62] and 
2,4-D [63]. The results of ODs calculation are summarized in 
Figure S3. For the studied compounds, the distances between the 
negative charges on the carboxyl groups and the partial positive 
charges on the rings ranged from 5.026 to 5.158 A and the 
maximal absolute values of point charges which are used for 
estimation of ODs are observed in HEA-pABA. 

Docking and QM/MM calculations 

Auxin's receptor, TIR1, belongs to a large family of F-box 
proteins, which function as part of a protein destruction machinery 
in the ubiquitin-proteasome system. At the bottom of the TIR1 
surface pocket, auxin helps nucleate a hydrophobic core together 
with TIR1 and its substrate polypeptide [31]. It is known that 
auxin molecules IAA, 2,4-D, and 1-NAA bind to the TIR1 pocket 
in a similar manner. Their common carboxyl group interacts with 
the same positive charged residue at the bottom of the TIR1 
pocket [55]. Therefore, in order to estimate the ability of the 
anionic form of p ABA (p AB A ) and pABA to interact with positive 
charged residue at the bottom of the TIR1 pocket, we performed 
the calculation of molecular electrostatic potential map on 
accessible surface area (Figure 3). In pABA the region of 
maximal attractive potential is located in vicinity of both oxygen 
atoms, while in pABA it is concentrated in vicinity of carbonyl 
oxygen. Thus, both oxygen atoms inpABA - are able to bind with 
positive charged residues at the bottom of the TIR1 pocket. In 
both molecules the attractive potential, which is located in vicinity 
of the rings, enhances their association with TIR1 pocket. 

The possible binding modes of p ABA to the TIR1 protein were 
then evaluated by performing docking experiments using the 
AutoDock 4.2 program. Both the neutral and the anionic form of 
p ABA were docked in the auxin binding site of TIR 1 protein as 
well as auxin itself (IAA), which was redocked to evaluate the 
performance of the program on this system. The best fit docking 
poses obtained from the simulation (see Materials and Methods) 
were then submitted to cluster analysis on a geometrical basis 
using a 2 A rmsd threshold. This clustering procedure revealed 
that in each case only one cluster is present, suggesting a very good 
convergence of the algorithm. In Figure 4, the docked structure of 
both neutral and anionic forms of pABA and of auxin IAA are 
shown. Protonation does not seem to affect binding geometry of 
pABA. It makes hydrogen bond with Arg403, Ser438 (using the 
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Figure 2. Crystal packing of HEA-pABA. A. Crystal packing of HEA-/ 
molecular packing in the crystal of HEA-pABA. 
doi:1 0.1 371 /journal.pone.01 01 892.g002 

carboxyl function) and with Leu439 backbone oxygen using the 
amino function. These interactions are the same polar interactions 
exhibited by IAA, even if in auxin the hydrogen bond of the 
indolic nitrogen is slightly shorter than the corresponding 
hydrogen bond in /ABA. 

Binding energy calculated by AutoDock for the best fit 
conformations are as follows (in kcal/mol): IAA, —7.91; pABA~ , 
— 6.19; pAHA, —5.97. Again, no significant difference exists 
between the neutral and the anionic form of pABA, while IAA 
exhibits a slightly stronger binding to the substrate, probably due 
to the larger size of the hydrophobic moiety. However, it is to be 
evidenced the quite small difference in binding energy provided by 
the three compounds, indicating a similar mode of binding of the 
proposed ligand with respect to auxin. Molecular dynamics 
simulations, on the other hand, showed that, while auxin and 
anionic form of p ABA conserve their position and orientation in 
the binding site compared to the docked structures, the neutral 
pABA is expelled from the binding site during the simulation. 
Accordingly, the corresponding Root Mean Square Deviations 
(RMSD) of the three ligands during simulation with respect to the 
MD starting structure show different trends (see Figure S4). While 
auxin IAA and anionic form of ^ABA display small and well 
equilibrated RMSD values, for neutral pABA the RMSD increase 
progressively up to more than 8 A after 10000 ps, indicating a loss 
of the starting structure as a consequence of its expulsion from the 
binding site. The MD results therefore suggest that the salt bridge 
interaction between Arg403 and the charged carboxyl group on 
auxin or pABA~ is crucial for the interaction of ligands with TIR1 
protein. 

Geometry optimizations at the hybrid QM/MM level were 
performed on the final MD structures of auxin and pABA 
complex. The QM/MM optimized structure of the TIR1 -auxin 
complex revealed no substantial modification with respect to the 
docked structures. The HB interactions between Arg403 and 
Ser438 and auxin carboxyl group are preserved, as well as the HB 
between the indolic nitrogen and the backbone carbonyl of 
Leu439 (see Figure S4). In the QM/MM optimised structure of 
the TIR 1 -p ABA complex, interaction of the carboxyl group with 
Arg403 and Ser438 is preserved, while no HB is formed between 
the amino group of f^ABA and Leu439 (see Figure S5). The 
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representing the 2-D layers parallel to (100) plane. B. Fragment of 




Figure 3. Mapping of the molecular electrostatic potential onto 
p-aminobenzoate anion (pABA - ) and pABA accessible surface 
area. 
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Figure 4. Binding modes of IAA (red), pABA (orange) and pABA" 

polar hydrogen atoms are shown. 
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(blue). Hydrogen bonds formed by pABA are drawn as green lines. Only 



binding energies evaluated at the QM/MM level are — 
1 17.0 kcal/mol and -85.6 kcal/mol for IAA and pABA - , 
respectively, in agreement with the observed differences in the 
network of interactions between the protein and the two ligands. 

It is worth noting that the larger magnitude of binding energies 
with respect to the values obtained by docking calculations should 
be ascribed to the fact that QM/MM calculations have been 
performed in vacuo, owing to the high complexity of the 
investigated systems. Within such approximation, the electrostatic 
interaction necessarily results to be overestimated with respect to 
the more realistic values obtained by docking experiments, where 
the solvent effect is implicitly taken into account by the score 
function used. In spite of such approximation, QM/MM 
calculations appear to confirm the similar interaction of anionic 
form of pABA with respect to auxin IAA. It is to be noted that very 
similar results have been recently obtained on the binding 
mechanisms of neutral and charged phenolic compounds with 
the human transthyretin, hTTR, using a similar approach 
combining molecular docking with QM/MM optimizations [64]. 
In agreement with their conclusions, ionization effects should not 
be neglected when modeling the binding modes of ligands with 
proteins. 

HEA-pABA and pABA activity on root growth 

Phytotoxicity tests, especially the seed germination and root 
elongation tests, have more advantages over those toxicity tests 
using animals and algae, such as sensitivity, simplicity and low cost 
[65-67]. Phytotoxicity of HEA-pABA, pABA and HEA were 
determined at five different concentrations. Seeds of the model 
plant Arabidopsis thaliana have been used for germination, early 
seedling growth and gravity response tests. Figure 5 shows wild- 
type (WT) A. thaliana root growth on a medium supplemented 



with 0.02, 0.1, 0.2, 1 and 2 mM of HEA-pABA, pABA and HEA 
as compared to controls after 5 and 10 days (see also Figure S6). 
HEA-pABA and pABA showed a clear difference in primary root 
growth, and lateral root development after prolonged growth (10 
days) (Figure 5). Moreover, the root growth difference was already 
apparent after 5 days, indicating that the action of tested 
compounds occurs throughout the period of elongation growth 
(see Figure S6). After 10 days, HEA treated seedlings exhibited a 
comparable effect to control for primary root length and an 
increased number of lateral roots at all concentrations tested. At 
the lowest concentration (0.02 mM), pABA treatment showed an 
increased number of lateral roots and primary root length, a 
common anti- auxin feature (Figure 5). At this concentration, 
HEA-pABA displays a comparable effect to control, regarding the 
primary root length and an increased number of lateral roots. At 
0.1 and 0.2 mM, both HEA-pABA and pABA promote lateral 
roots, while they inhibit primary root length (Figure 5). The 
phenotypes like short root length, increased lateral root number 
are phenotypes similar to those observed in A. thaliana treated 
with auxins [68-70]. At high concentration (2 mM), primary root 
length is inhibited but also displayed a different growth behavior 
with respect to control seedlings. Preliminary tests aimed to 
compare HEA-pABA, pABA and IAA effects indicate that the 
latter exerts a 1000-fold higher root length inhibition (data not 
shown). The agravitropic effect of root growth is specific forpABA 
and HEA-pABA, never having been reported for an auxinic 
compound. Therefore, we argue that HEA-pABA and pABA 
might potentially act as auxin-like regulators of root development 
in A. thaliana (0.1 and 0.2 mM) and display an agravitropic root 
response at high concentration (2 mM) (see Figure S7). 
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Figure 5. Root length and number of lateral roots after 10 days of treated A. thaliana seedlings. HEA-pABA, pABA and HEA treatments at 
different concentrations in comparison with control are shown. Control corresponding to pABA was noted CONTROL 0.2% EtOH and for HEA-pABA 
and HEA, CONTROL. Values with different letter annotation are significantly different (P<0.05). Data are means ± SE of 5 replicates. 
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Conclusions 

This study presents the synthesis and the unique structural and 
biological properties of HEA-^)ABA. In the crystal, the salt unit of 
HEA-^ABA serves as building block of the supramolecular 
architecture and in the crystal it is self-assembled via ionic N- 
H " O and normal O-H "O hydrogen bonds. The pyramidal 
geometry of jf)ABA amino group allows to form three hydrogen 
bonds in the crystal. The transfer of proton in HEA-^)ABA has led 
to significant polarization of anion in comparison with jbABA and 
the absolute values of point charges, which are used for estimation 
of ODs, are higher in HEA-^)ABA. For both these compounds, 
ODs are equal to 5.045 and 5.026 A respectively. Both oxygen 
atoms in anionic form of /)ABA are able to bind with positive 
charged residues at the bottom of the TIR1 pocket. In //ABA and 



pABA the molecular electrostatic attractive potential, which is 
located in vicinity of the rings, enhances their association with 
TIR 1 pocket. Docking experiments on complexes of p AB A (in 
both neutral and anionic form) and IAA with TIR1 provide a 
similar binding mode between the two ligands. MD and QM/MM 
calculations confirm such result as far as jbABA - is concerned. 
Phytotoxicity test indicates that compound HEA-/) ABA and jbABA 
act as auxin-like regulators of root development in A. thaliana (0.1 
and 0.2 mM) and display an agravitropic root response at higher 
concentration (2 mM). This study suggests for these compounds a 
potential plant growth regulatory activity of the auxin type. 
Between the two ionic species, anionic from acid and cationic from 
HEA, the anionic species presented an increased phytotoxic effect. 
However, up to now, the plant growth regulatory properties of 
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HEA-^ABA and/)ABA have not been found in literature, which is 
needed in order to develop its new application fields and carry out 
relevant theoretical research. 

Supporting Information 

Figure SI UV absorption spectra of HEA-pABA and 
pABA in 0.1 M NaOH solution. 

(PDF) 

Figure S2 FT-IR spectra of HEA-^ABA (black) and pABA 
(red). 

(PDF) 

Figure S3 Optimal distance calculation for HEA-/»ABA, 
pABA and classical auxin molecules (IAA, 1-NAA, 2,4-D). 

(PDF) 

Figure S4 (A) Binding mode of IAA with TTR1. Hydrogen 
bonds are drawn as green lines. (B). Root Mean Square 
Deviation of the ligand with respect to the MD starting 
structure. 

(PDF) 
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